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Abstract 

We study disk amplitudes whose boundary conditions on matter configurations 
are not restricted to homogeneous ones. They are examined in the two-matrix 
model as well as in the three-matrix model for the case of the tricritical Ising model. 
Comparing these amplitudes, we demonstrate relations between degrees of freedom 
of matter states in the two models. We also show that they have a simple geometrical 
interpretation in terms of interactions of the boundaries. It plays an important role 
that two parts of a boundary with different matter states stick each other. We also 
find two closed sets of Schwinger-Dyson equations which determine disk amplitudes 
in the three-matrix model. 
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1 Introduction 


It is well known that (m, m + 1) nnitary conformal model can be constrncted microscopi¬ 
cally as a lattice statistical model[|^. At each site this statistical model has local degrees 
of freedom labeled by the points of A^-i Dynkin diagram. The (m, m-|- 1) model coupled 
to 2d gravity can be described by matrix model. The {m,m -|- 1) model has {m — 1) 
microscopic degrees of freedom and the (m — l)-matrix chain model|^ naturally corre¬ 
sponds to the (m, m -|- 1) model coupled to gravity. On the other hand the two-matrix 
model P, ^ H can also describe this system near an appropriate critical point though 
there are only two matrices. In this paper we address the correspondence between the 
degrees of freedom described by the matrices in the (m — l)-matrix chain model and those 
in the two-matrix model[0. 

As a non-trivial simplest case, we study the (4, 5) tricritical Ising model coupled to 2d 
gravity. As actions of the two- and the three-matrix models which realize this system, we 
take 

S{A,B) = Etr{u(A) + U{B)-AB}, (1.1) 

S{A,B,C) = jtr {Ui(.A) + UiiB) + Ui(.C) - AB - BC] . (1.2) 

An amplitude of loops which have homogeneous matter conhguration corresponds to 
an expectation value of fr(A") and so on. The integration of the matrices is to be separated 
into one over eigenvalues and one over angular variables. We can integrate the angular 
variables hrst and reduce the original integral to that in terms of the eigenvalues. Then 
through the orthogonal polynomial method we can show that the loop composed of the 
matrix A{B) in the two-matrix model corresponds to that composed of the matrix A{C) 
in the three-matrix model. On the other hand, an amplitude of loops with heterogeneous 
matter conhgurations corresponds to an expectation value of tr{A"'B^ ■ ■ ■) and so on. We 
cannot integrate the angular variables hrst in this case. Then the argument in the case 
of homogeneous loops can not be applied, and the correspondence between the matrices 
in the two models is not so trivial in this case. For simplicity we restrict our attention to 
disk amplitudes. One of the purposes of the paper is to calculate disk amplitudes with 
heterogeneous loops in the two models using Schwinger-Dyson technique^, ||, ^ |TI[, and 
to study the correspondence between the matrices in the two- and three-matrix models. 

Studying these disk amplitudes is also a very interesting problem by itself. We are 
forced to study heterogeneous loop when merging different homogeneous ones into a single 
one. The amplitudes involving heterogeneous loops have not been studied deeply^]. 

^ In the case of Ising model these amplitudes were studied in §0- 
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We will obtain an interesting geometrical pictnre on the amplitndes involving hetero- 
geneons loops in sect.4. If two parts of a loop have different matter states, these are forced 
to stick each other and the original heterogeneous loop reduces to homogeneous ones. We 
obtain this picture through the study of the disk amplitudes. We believe that this picture 
can naturally be extended to cases with many loops or cases with many handles because 
the sticking of boundaries of loops is a local phenomenon and does not depend on the 
global nature of surfaces. 

The paper is organized as follows. In sect.2 we compute heterogeneous disk amplitudes 
whose loops are composed of two arcs using the Schwinger-Dyson technique in the two- 
matrix model. In sect.3 we compute similar amplitudes in the three-matrix model. In the 
process of our calculation we hnd two closed sets of Schwinger-Dyson equations which are 
composed of seven and ten equations respectively. The success of our calculation is due 
to these hndings. In sect.4 we discuss the correspondence between the matrices in the 
two- and the three-matrix models, and provide a geometrical picture on these amplitudes. 
Sect.5 includes a summary and discussion. 

2 Two-matrix model case 

As a critical potential which realizes the (4, 5) model in the two-matrix model, we take 

U{(f)) = 80-|-20^-|--0^-|--0^, (2.1) 

which can be determined by the method of orthogonal polynomial [^] (see appendix). In 
this section, we would like to calculate the disk amplitude 

All “A 

= >= £ V”-'. (2.2) 

N p- Aq- tS N 

and its continuum universal part WAB{CA,CB,t) in the large N limit by means of the 
Schwinger-Dyson technique. The boundary is consisted of two parts which have different 
matter states. 

Let us work out the disk amplitude 

OO A 

w'f(p) = E 77 (2.3) 

n=0 

which will be necessary for the calculation of the amplitude ( p.2D . Consider the Schwinger- 
Dyson equations 

0 = Y. J[dAdB]-^ {tr , (2.4) 
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where we decomposed the matrices as A = Yla=i ^tc. by introducing the basis {t“} of 
the hermitian matrix. Using a notation = -^{tr{A'^B^)), we can rewrite eqs. (|2.4| ) 

as 

n—1 

0 = - 8[A^B^] - - 8[A^^^B'^] - [A^+^B'^] + 

1=0 

0 = 8[A”] + 4[A’"S] + 8[A^B‘^] + [A'^B^] - [A”+^], (2.5) 


in the large N limit. It is convenient to use the resolvent representation eq.( p.3|) and we 


obtain 


0 = {WA{p)-x{p)}wf\p) + wf^"\p) + a^^"\p), (2.6) 

0 = (S-p)WA{p)+iW[^\p) + 8Wi^^'\p) + w[^'\p)+A, (2.7) 

where 

x{p) = 8 + 4p + 8p^+p^, 

a^^'\p) = {4: + 8p + p^)[B’^] + {8+p)[AB'^] + [A^B’^]. (2.8) 

Note that [1] = A and we used the Z 2 symmetry. One can easily hnd that eqs.(^.6|) for 
fc = 0,1,2 and eq. (p?7|) make a closed set of equations |^, ^]. We can eliminate lU 4 ^^(p), 
IU 4 {p) and IU 4 (p), and have the following fourth order equation of Wa{p)- 

VAip)"^ + a^{,p)VA{,pf + a2{p)VA{pf + ai(p)14(p) + ao(p) = 0, (2.9) 

where 


Va{p) = Wa{p)-x{p), 

« 3 (p) = x{p) - 8, 

a 2 {p) = A - 8x{p) + a^^\p), (2.10) 

«i(p) = p — 8 + 4a;(p) — 8a*'^^(p) — a*''®^(p), 

<^o(p) = — A + (p — 8)a;(p) + 4a*^^^(p) + 8a‘''®^(p) + a*''®^^(p). 

We must provide the amplitudes [A], [A^], [AH], [A^H] and [A^H^j in order to solve 
eq. &■ These can be determined by the method of orthogonal polynomial (see ap¬ 
pendix). The continuum limit can be carried out by the renormalizaion A = 70 — lOa^t 
and p = 2a(A with the lattice spacing a [^. Here 70 is a critical value of A. Assuming 
the scaling form of U 4 (p) as 

Va{p) =Co + CiCAa + C2 Wa{Ca, + C>(a®/^), (2.11) 
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and substituting this form into eq. (^^) , we have the equation of WA{CA,t), 

wa{Ca, t)^ - ^t^^^WAiCA, tf + ^ 20fCl - 16Cl) = 0 , (2.12) 


and find (co,Ci,C 2 ) = (0,-2, ±2). By solving eq. (|2.12| ), we can obtain the continuum 
universal disk amplitude WA{CA,t) as 


WA{CA,t) — + \lC a~'^ + ~ \JCa~'^ 


(2.13) 


We examine next the amplitude Wab{Pi q) which is the prime interest. Observing that 
^ab{p, q) = \p)(l~’^~^^ from eqs. (p. 6 |) , ( p.7|) we can find the following equation, 


^ (4 + 8 p + p^-)WA{q) + (8 + p)w!P{g) + wf\q) - Wa(p) „ 

’ x{p) -q- Wa{p) 


A careful consideration is needed to extract the universal amplitude WAB{CA,CB,t) from 
this equation. For example let us consider (p) for finite m. The boundary of the 

corresponding disk involves a part of finite lattice length composed of the matrix B. Then 
the contribution from such a part in Wab{p, q) turns out to be non-universal. In general, 
any polynomials of and (b multiplied by IF 4 [p) or hFg ’{q) are non-universal 
quantities of 1 F 4 b(p, <?). Polynomials of (a and (b are also non-universal. We should 
drop these quantities appropriately to extract a universal part of Wab{p,q)- Using the 
expansion of Va{p) (| 2 . 11 D and a similar expansion of VB^q), we can find 


Wab{p, q) - 2(4 + aCA + <b) {Wa{p) + WB^q)) 

- (F®’(p) + wTl?)) 

= “7 — 32(Ca -|- C,b)0‘ — 4(Cl + 4 CaCb + 

+4 -t- wa{C,a, t)wB{C,B, t)) . (2.15) 

In the left hand side, we appropriately subtracted some non-universal quantities in ad¬ 
vance. Moreover we should drop any terms which are analytic in both (a and (b from 
the right hand side. Therefore we can read the continuum universal part of Wab(p, q) as 




( 2 . 16 ) 


where WBiC^t) = 'a^A(CU) from the Z 2 symmetry. It should be noted that the terms with 
order higher than Va{p) ( | 2 . 11 |) do not appear in the right hand side of eq. (p.l 6 ), 

so that wab{Ca, Cb, t) can be expressed only in terms of wa{Ca, t) and wb{Cb, t). We will 
discuss the implication of this fact in sect.4. 
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3 Three-matrix model case 


In this section, we will investigate the disk amplitudes 


_ “A 

»'AB{p,g)= T. 

n,m=0 


A 


—n—l^—m—1 


WAc(p,r)= E ^{«r(yl”C”'))p-”-V 


n,771=0 


(3.1) 

(3.2) 


and their continuum universal parts WAB{CA,CB,t), WAciCAXcyt) in the three-matrix 
model. As potentials which describe the (4, 5) model, we take 


[/i(0) = 


111 

Iq' 

3 




Vs 


12 


(3,3) 


These can be found by the orthogonal polynomial method (see appendix). In order to ob¬ 
tain the amplitudes (^T]), (|0[) , we have to calculate ^ \p) = 
and where [A^B^C^] = ^{tr{A^B^C^)). 


Let us examine W\ (p) first. Consider the Schwinger-Dyson equations 

d 


0 = j[dAdBdC]^[tr(A^eB^C^)e-^^^’^^^'>], 

0 = E/ [d^dBdC]^ {tr (A^t^C^) , 


(3.4) 


0 = E/ [dAdBdC]^ [tr {A^B^t'^) e-5(^-S,c)| _ 

We may write them in the resolvent representation, 

0 = {ir^(p)-9(p)}w_f”°*’(p)+'ir!f"‘"'‘’(p) + a'®"^'‘’(p), 

0 = -|irW(p) - - wf^lp) - pwrip) + iA% 


0 = ■jj'^ P (P) - 2 ^^ A (P) - (P) - ^'^A (P). 


(3.5) 

(3.6) 

(3.7) 


where 


111 9 2 

y{p) = “ 2 ^ ~ ^ ’ 

Here [1] = A and we used the Z 2 symmetry. One can hnd that eqs. (El) for (m, k) = 
(0, 0), (0,1), (1, 0), (1,1), eqs. (P^) for fc = 0,1 and eq. ([H7^) for m = 0 are independent and 
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make a closed set of equations^]. By eliminating \p), 

(p) and (p), we obtain the following fifth order equation of WAip), 

UA{pf + ai{p)UA{pY + ^sipWAipf + a2{p)UA{pf + ai{p)UA{p) + ao(p) = 0. (3.8) 

Here 


Ua{p) = WA{p)-y{p)^ 

a^{p) = -12 + p(p), 

a^{p) = 18 + 8p + a^^\p) - 12?/(p), 

«2(p) = 92 — 48p — 12a^^^(p) — a^^^(p) — 4A + (18 + 8p)p(p), 

«i(p) = —111 — 72p + 16p^ + (18 + 4p)a‘'^^(p) + 6a^^^(p) — 4a^'^^(p) 

+24A+(92 -48p)p(p), (3.9) 

«o(p) = (92 — 24p)a*^^^(p) + (18 — 4p)a^^^(p) + 24a*^'^^(p) + 4a^^‘^^(p) 

+ (72 - 16p)A + 16[A] + (-111 - 72p + 16p^)p(p). 

The amplitudes [A], [H], [AH], [AC] and [AHC] are determined by the orthogonal poly¬ 
nomial method (see appendix). With the renormalization, A = 35 — |a^t and p = \a(^Ai 
we assume the scaling behavior of Ua{p) as 

Ua{p) = Co + ciC^a + C2Wa{Ca, + 0(0®/"^). (3.10) 

Substituting eq. (|3.10|) into eq. (|3.8|) and after similar calculation in sect.2, we hnd that 
(co, Cl, C2) = (0, 2, ±2/3) and 

WAiCAit) = + \jCA — \/Ca~ ■ (3-11) 

As expected, this coincides with the result for wa{Ca, t) in the two-matrix model, eq. (^.13D . 

Next let us examine Wb{(i)- In this case, we found that ten Schwinger-Dyson equations 
are needed. For example let us consider the following ten equations: 

0 = ^ J[dAdBdC]-^{tr{t^B^)e-^^^’^’^^} , 

0 = ^y'[dAdHdC]^{tr(ACH")e-^(^’^’^)}, 

0 = ^J[dAdBdC]-^{tr{Ct^AB^)e-^^^’^’^'>}, 

0 = ^y'[dAdHdC]^{tr(CACH’*)e-^(^’^’^)}, 
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0 


(3.12) 


= E / [dAdBdC] ^ {tr (ABCeB^) } , 

0 = Y.J[dAdBdC]^{trieB-)e-A^^^’^)]^ 

0 = ^J[dAdBdC]-^{tr{{CAe-Ct‘^A)B^)e-^^^’^'^'>], 

0 = ^J[dAdBdC]-^{tr{CBBB^)e-AAB,c)y 

0 = ^J[dAdBdC]-^{tr{{ACAt'^-ACBA)B^)e-A^'^’^^], 

0 = ^J[dAdBdC]-^{tr{{CBAt‘^-CBeA)B^)e-A^’^’^)y 

In the resolvent representation, we have 

0 = {WB(?)-2(?)}irB(<;) + 2TF*)'’(<;) + i)''>(g), 

0 = {WB(rD-4q)}w^B\<i) + ^A('D + w'A\'t} + i^^\'D, 

0 = {WB(q)-z{g)}W'B''\q)+'i^f'"\q) + w'B'^^\g) + b'--*'^\,]), 

0 = {WB{q)-z(q)}w‘B‘’°\q) + W^B'^^*\q) + w'B''^-\q) 

+|3l|Tri5’'(,)+()<'’®CI(9). (3.13) 

0 = {^-q)wB{q)-lw^B'(q)-wA’(.q)+A, 

0 = qW''i^'{q)-W'B‘‘'^'’{q)-[ACl 

0 = {^q-q^)w'f{q)-lw‘^A''hq)-w'f‘’^''(.q)-(^-q)lA] + lAB], 

0 = \AC]WB(q)-\A]w^B\q)-q^A^\q)+'^f'^\q) + \A"c], 

0 = qW'f°\q)-WT'‘'\q)-[ABC], 

where 

/ \ 3 1 2 

A<i) = 

= -{^ + ^q)[A^B^C^]-^[BA^B^C’‘] . (3.14) 

We can hnd that these are independent and make a closed set of equations. The fourth 
order equation which determines Wsiq) is obtained as: 


UB(qY + Mq)UB(qf + Mq) = o. 


( 3 . 15 ) 



where 


Usiq) = WB{q)-z{q)- 


9 


The coefficients of this equation are given by 

A 


/^2(?) 

^ 0 ( 9 ) 


.HI2-3A-H 

4 2 


/35 

+ (y“ 2 


33 1 4 

^- 4 ^ -16^’ 


972 + - 46[A] -%[B] + ^A[B] - 6[AC] + - 9[A5] - [ABC] 

8 4 16 4 16 

+ {-729 + - [AC] + ^[B] + iA[B] + |[AB]} , 

[4B] 1 , 
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[^] 


qo 


-T 


(3.16) 


With the renormalization q = 2a(B, solving eq. (|3A5|) directly, we can hnd the the disk 
amplitude UB^q) as 

UBiq) = ± WB{CB,t)a^^^ + 0{a^/^) , (3A7) 

where WBiCsyt) coincides with WA{CB,t)- 

Now let us turn to the calculation of W ab{p, q), W Acipc)- From eq. fP.Sl) for k = 0, 
one can obtain the relation 


, > (l+PWB{q)+w‘f{q)+WA(p) 

WAB{p,q) = - -=----y—-. 

Wa{p) - y[p) + q 


(3.18) 


By combining eq. (|3A8| ) and the hrst equation of ( PA3| ), WAB(yP,q) can be expressed 
in terms of Wa{p) and WB^q)- In order to extract a universal part, we must drop 
polynomials of (a and (b multiplied by Wa (p) or Wb (q) as well as polynomials 
of both (a and (b appropriately, because of the same reason as stated in sect.2. Using 
the expressions ( p.lOD and (|3.17|) , we can hnd 

WABip,q) - ^IFa(p) - ^WBiq) 


1607 9 , , 

^ + ^(llCA + 8CB)a 


(3.19) 


1wa(CaA)^ + 2m;a(CaA)ws(CsA)+ 2M;s(Cs,t)^-2CA ^ 

8 Ca + Cb “ + 

In the left hand side, we subtracted some non-universal quantities in advance appropri¬ 
ately. Moreover we should drop hrst and second terms in the right hand side, because 
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they are polynomials of both (a and (b- From this equation, we can hnd the continuum 
universal disk amplitude WAB{CA,CB,t) as 


wab{Ca, Cb, t) 


u!a{Ca, + 2wa{Ca, t)wB{CB, t) + 2wb(Cb, 

Ca + Cb 


(3.20) 


We can observe that terms with order higher than in Wa{p) and Wb{(i) do not 
appear in the right hand side of eq.( p.20|) . Thus wab{Ca, Cb, t) is expressed only in terms 
of wa{Ca, t) and wb{Cb, t). 

Next let us consider WAcip,'’^)- From eqs. (|3.5|) for m = 0,1 and eq.(|3.6|), we can 
obtain the equations 


WAc{.P,r) 


(I + p)Wc{r)_+ W^c\r) - trig(p,r) 

'/Va(p) - v(p) 


(3.21) 




(I +p)wTV2_+ tTc'’'(0 - W^Ac(P,r) 

Wa{,p) - y{,p) 


(3.22) 


\w^Ac{P.r) + Iw^A^iP^r) + ip + r)WAc{p,r) - Wa{p) - Wc{r) = 0 , (3.23) 


respectively. By combining these, we can express WAc{p,r) in terms of Wa{p) and 
Wcij"). Note that Weir) = W a{j) because of the Z 2 symmetry. Using the expression 
( p.lOj) and (|3.17|) , we hnd 


41 


ic^). 




WAcip,r) - ^ hF)r \p) + ^(r) - 


801 

112 


W‘-P(p) + W'c'{r 


ic), 


aCc) Wa{p) - (^y - Wc{r) 

1 


90 _ 11 

y ~ ^ 

44649 99,, , . 

’ 7168 + 

73 


90 11 


3584 


(-23200f + 1377C1 + 8120 CaCc + 1377C 


+ + q'^a{.Ca, t)wc{Cc, t)j + C>(a^^/^) . 


(3.24) 


Here we subtracted some non-universal quantities in advance appropriately from WAcip, r) 
The hrst, second, third and terms in the right hand side should be dropped, because 
they are polynomials of both y and Cb- We can read, therefore, the continuum universal 
disk amplitude WAc{CA,Cc,t) as 


wac{Ca, Cc, t) = wa{Ca, t)wc{Cc, t) . (3.25) 
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4 Comparison and interpretation 


In the previous two sections, we obtained the disk amplitudes with heterogeneous bound¬ 
aries WAB{CAXB,t), wabXaXbX) and WAcXAXc.t) ( eqs. (|2.iq) , (|3.2q) and (|3.25|) re¬ 
spectively). In this section, we will compare them and provide a geometrical interpretation 
of these amplitudes. In this and the next sections, we denote the matrices A, B and C 
in the three matrix model as A, B and C respectively, in order to distinguish from those 
in the two-matrix model. We will refer to a part of boundary which is composed of the 
matrix A as “boundary ^4” and so on. 

From eqs. (|2.16|) and (|3.25|) , we observe that wab and wac have exactly the same 
form. We can consider that boundaries A and B correspond to boundaries A and C 
respectively. In the case of loops with homogeneous matter states, this correspondence 
is natural from the view point of the orthogonal polynomial method. In the case of 
heterogeneous boundaries, however, the method of the orthogonal polynomial cannot be 
applied and this correspondence is not so trivial. 

From eqs. (|3.20|) and (|3.25| ), we can observe that wab and wac have quite different 
forms. But wab obtained by eq.( |3.2U|) has a very similar form to the disk amplitude 
in the case of the Ising model realized by the two-matrix model 


, 11 


WabXaXbA) = 
Here is given by 


w^^\CA,ty+w^^KCA,tW^KCB,t)+w^^KCB,ty - 

Ca + Cb 


(4.1) 


/ /-\ 4/3 / , -\ 4/3 

w"' (C, i) = (c + /c" - 1 ) + (c - Vc" -«) 


(4.2) 


Note that, however, there is no symmetry under interchange of A and B in eq. ( |3.2(j| ). 

For the sake of discussing why wab and wac have so different forms, it is useful to 
consider the inverse Laplace transformed amplitudes. Let us denote the inverse Laplace 
transformed amplitudes of WA{CA,t) as Wa^X) etc. For example, wabXaXbA) and 
^ab{(^aXb) are related by the equation 


poo poo 

'WAB{CAXB,t)=CA^B[y^AB{^AXB)]= Ma f's)- (4.3) 

JO JO 

Here WABi^A^B) represents a disk amplitude where length of each part of the boundary 
is fixed. First we easily obtain the relations: 


mB^A^B) = Wa(^a)Wb(^b) , (4.4) 

^^AciXXc) = W^(£^)Wc(£c) • (4.5) 
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As for we use the following formulas of the inverse Laplace transformation; 


L 


-1 


= e 


—at 


Qj 

£-1 [e-^“F(C)] = B[l - a) f{i - a) , 

'F{cy 


C 


-1 


c + ® 


= e-“^ / dxe^^fix) , 
Jo 


(4.6) 

(4.7) 

(4.8) 


where F{() denotes the image of the Laplace transformation of f{i)- Using the formulas 
( ^(^ and (|4)7|) , we obtain the relation: 




v^a((aA)^ 

= [ 

Ca + Cb 

A [ 


= 0{iA-^B) F/ IwAiCA^t) 


= 0{iA - 4) 


A 

^A~^b 


/o 


dmA{^)y^A{^A-^B-^) 


We also have 


C~FC~F 


Wa{Ca^ J)wb{Cb-: t) 


Ca + Cb 


= C~F 


die-M‘A-e>WAi)wB(<;B,t) 


re A 


diWAimis - 4 + ^)>Vb(4 - ^4 + ^) , 

JO 

by using the formulas (^4.8|) and ([4.7|) . From (^4.6|) we have 

rfC-A =C-A[e-‘^>‘A\=6{i^-iB) ■ 

Collecting eqs. (^^) - ([4.11|) together, we obtain the expression for Wab) 

y^ABi^A^B) = 0{^A - ^b) M WAimAi^A - 4 - ^) 

Jo 


+2 


+2 


di w^(£^-£)Ws(4-^) 


B 


-h) 


'eB~eA 


di >Vb(^)>V 5 ( 4 -^ 4 -^) 


- 4) . 


(4,9) 


(4.10) 


(4.11) 


(4.12) 


Now let us consider the geometrical meaning of eqs. (^4.5|) and ( |4.12| ). As for eg. ( |4.5|) , 
it is easy to understand that a loop composed of boundary A and C splits into two loops 
with different homogeneous matter states (see fig.|^). On the other hand, the hrst term 
in the right hand side of eg. (|4.12|) represents the conhguration depicted in £g.^(a). All 
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A 


C 





Figure 1: The original loop composed of two different parts of boundary splits into two 
loops each of which has homogeneous matter conhgurations. 




(a) (b) 

Figure 2: Due to the sticking of two different kinds of boundaries, the original loop splits 
into two loops with homogeneous matter conhgurations. 
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region of the boundary B is stuck to the boundary A, and the original loop also splits 
into two loops with homogeneous matter states. Likewise the second term in eg. (|4.12|) 
corresponds to the case in £g.^(b). Parts of boundaries A and B are stuck each other, so 
that the original loop splits into two loops with homogeneous matter states. The fourth 
term represents the contribution from the case where the boundaries A and B are stuck 
completely. From this geometrical picture, we can conclude that the original loop, in fact, 
splits into two loops with homogeneous matter states. 

Next we consider the reason why there is a difference between Was Wac from 
this geometrical point of view. In the case of Was fhe boundaries A and B stick each 
other. On the other hand, the boundaries A and C are connected at only one point in 
Wac- This difference can be accounted for as follows. The (4, 5) minimal conformal model 
has three degrees of matter freedom labeled by the points of the A^ Dynkin diagram. We 
can interpret that the matrices A and C correspond to the ends of the diagram and B to 
the middle point. The boundaries A and B stick each other, because the corresponding 
states of A and B interact directly as opposed to A and C. In the case of the two-matrix 
model eq.( ^4.4|) , the matrices A and B correspond to the ends of the Dynkin diagram. The 
boundaries A and B, therefore, do not stick each other. 


5 Summary 

In this paper, we have considered the (4, 5) minimal model coupled to 2d gravity described 
by both the two- and the three-matrix models. We have calculated the disk amplitudes 
with non-trivial boundary conditions for the matter conhgurations, and have shown ex¬ 
plicitly the relation among the matrices of these two models. A geometrical interpretation 
of the resulting amplitudes have also been obtained. In the process of our calculation in 
the three-matrix model, we found that seven and ten Schwinger-Dyson equations make 
two closed sets. These two sets of equations determine the disk amplitudes Wa{.p) and 
W§{q) respectively. 

We obtained the universal disk amplitudes wab, was and Wac^ whose boundaries 
composed of two arcs of hnite length with different matter states, as eqs.( p.l6D , (|3.2CI|) 
and (|3.25|) respectively. We learned that the matrices A and B in the two-matrix model 
correspond to the matrices A and C in the three-matrix model. 

The geometrical meaning of these results is that the loop of the disk W as or ^Ac splits 
into two loops with homogeneous matter states. Only these conhgurations contribute to 
the amplitudes. 

In this paper, we only studied disk amplitudes with two arcs. It is straightforward 
to generalize our calculation to disk amplitudes with more than two arcs. Using the 
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Schwinger-Dyson equations, we can compute recursively the disk amplitudes with n arcs 
from those with smaller numbers of n arcs. 

What extent can we apply the technique in this paper to other matrix models? The 
point is whether we can construct closed sets of equations which determine the disk 
amplitudes with homogeneous boundaries (ex. W and W^{q)). In the three-matrix 
model of closed chain type, we can hnd that this is the caseQ We believe that the 
Schwinger-Dyson technique could be applied successfully to the multi-matrix models of 
closed chain type as well as open ones. 

We believe that the geometrical picture in this paper can be extended to the cases of 2d 
surface with many loops and handles. This expectation is natural because the interaction 
of boundaries is a local phenomenon and does not depend on the global nature of surfaces. 
This picture must be applied also to the case of general unitary minimal matter. Let us 
consider the (n -|- 1, n -|- 2) unitary model coupled to 2d gravity realized by the n-matrix 
chain model. The matrices naturally correspond to the matter degrees of freedom labelled 
by the points of An Dynkin diagram. We expect that two parts of the loops stick each other 
if the two corresponding matter states are adjacent in the Dynkin diagram. Due to this 
mechanism, heterogeneous loops must reduce to homogeneous ones. This phenomenon 
may be related to the formulation of the non-critical string held theory |T3|, pT[] , which is 
constructed in the limited space of loops with simple matter conhgurations. 
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Appendix A Orthogonal Polynomial Method 

In this appendix, we derive the critical potentials (p.l| ), ( |3.3| ) and critical values of A by 

following ref. [^]. We also show how to evaluate -^{trA^) etc., which are necessary to solve 

the Schwinger-Dyson equations in the text. We show the details of the calculation by 

restricting our attention to the tricritical Ising model; {p,q) = (4,5). 

^For the action S = N/Atr{V (T) -|- V(B) + V(C) — AB — BC — CA} with third order potential V, we 
can find that eight Schwinger-Dyson equations make a closed set. These leads to a fourth order equation. 
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First, we consider the two-matrix model. The potential U{(p) in the action (|1 . 1|) is an 
arbitrary polynomial; 



k=i 

(A.l) 

The integral over matrix elements can be converted into the one over 

the eigenvalues: 

Z{A) = 

f Y[dxidyiA{x)A{y)e 

i=l 

(A.2) 

S{xi,yi) = 

N 

— {U{xi) + U{yi) - xiyi) , 


A{x) = 

- Xjf- 



i<j 


By introducing the orthogonal polynomials n„(a;) which satisfy 


{m\n) = 

J dxdye = Smn, 

(A.3) 

we denote the matrix elements 

as 


Aimn 

{m\x\n), Pmn = {m\—■ —\n}. 

N ox 

(A.4) 


We can also derive the eqnation of motion 

Pmn = {m\U'{x)\n) - Xnm- (A.5) 


In the large N limit, the matrices X and P are replaced with the classical fnnctions 
X{z,X) and P{z,A) respectively: 

P{z, A) = U'{X{z, A)) - XiXjz, A). (A.6) 

Let ns determine the critical potential which realizes the (4, 5) model. We know that 
X will be a fourth order differential operator in the continuum limit. At the critical 
point, therefore, we can set X{z^A = A,,) = (1 — zY'jz and P{z,A = Ac) = AcZ+ (higher 
powers of z). After substituting these into the equation of motion, we can hnd the critical 
potential 

=8 + 40 + 802 + 03, (A.7) 


and the critical value of the cosmological constant as Ac = 70. 

Off the critical point, we set the classical functions X{z,A) and P{z,A) as 

X{z,A) = ^R{A)/z + a(A) + b{A)z + c{A)z^ + d{A)z^, 

P{z, A) = Az/\/R{A) + {higher powers of z). 
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Here -R(A) is called the specific heat function. Expanding the equation of motion in 
powers of we obtain the third order equation of a(A): 

a(A)3 + A)a(A)2 + m{R, A)a(A) + n{R, A) = 0, (A.8) 


where 



i{R,A) 


m{R, A) 


n{R, A) 

functions 

are given in 

6(A) = 

■/m) 

1 - 3R{A) 

c(A) = 

R{A){S + [ 


3 + 271i?(A) + 9i?(A)2 - 21R{kf 
1 + 15i?(A) ’ 

8 + 40i?(A) + 24i?(A)2 - 72R{Af 
1 + 15i?(A) ■ 


A = 


4)}, d{A)=R{A)l, 

3R{Af + {l8a(A)2 + 96a(A) + 128} R{Af + {3b{Af - l} R{A) 
+ {3a(A)2 + 16a(A) + 4} b{A)^R{A). 


The third order equation (|A.8|) has three possible solutions. Generally, the {p, q) model 
has a relation between the cosmological constant A and the specific heat R{A): 


R{A) - 1 ~ (A-A,)2/(p+9-i). 


(A.9) 


We should take the solution that satisfies this relation for the case of the {p, q) = (4, 5) 
model. 

The exact expression of X{z,A) can determine the expectation value -^{trA^). In the 
large N limit, the summation is replaced with the integration: 

A , A , A r dz 




(A.IO) 


Next, let us consider the three-matrix model. As in the case of the two-matrix model, 
we introduce the orthogonal polynomials n(a;) which satisfy 


(m|n) = / dxdydze 


-S{x,y,z) 


n(2:)n(a:) = d^n- 


(A.ll) 
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It is useful to introduce matrices Xi, X 2 , Xs and Pi; 



{m\x\n), 


{m\y\n) 


{m\z\n) 

[-Pllmn = 

, , A d 

<”liv ■ S 


Inm ’ 


With these matrices, the equations of motion are expressed as 


|U1„„ = {MU[(x)\n)-\X^]^^, 
(m|l'2fe)|n) = l-'-'llmn + l-'-'llmn- 


(A.12) 


Introducing classical functions, in the large N limit, (|A.12D can be rewritten as 


P(^,A) = P{(Xi(;^,A))-X 2(;^,A), 
P'(X2(^,A)) = Xi(z,A)+Xi(l/z,A). 


(A.13) 

(A.14) 


Remark that the classical function satishes X 2 {z, A) = X2{1/z, A), because the matrix X2 
is symmetric under the transposition. Now let us determine the critical potentials which 
realize the (4, 5) model. The classical functions are now 


Xi(z,A,) 


X2(^,AJ 


{l-z/m-zY 
(1 - zY 


By using these critical behaviors and the equations of motion (|A.13|) , ( |A.14|) , we obtain 
the critical potentials and the critical value of the cosmological constant. 


U'M) = 

Ac = 35. (A.15) 


The critical potentials determine the classical functions Xi{z,A.). We expand them in 
terms of 

Xi(^, A) = R{X)/z + a(A) + h{X)z + c{X)z^ + d{A)z^ + ^(A)^"^, 

X 2 {z,A) = f{A)/z‘^+ g{A)/z + h{A)+g{A)z +f{A)z‘^, 

P{z,A) = Az/\JR{A) + {higher powers of z). 
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After substituting these into (|A.13|) and (|A.14|) , we get the fourth order equation of a(A): 
a(A)^ + k{R, A)a(A)^ + i{R, A)a(A)^ + 7h{R, A)a(A) + n{R, A) = 0, (A.16) 


where 


k{R,A) = 

288R{Af + 60R{Af + 9 

32i?(A)3 + 1 

i{R,A) = 

12R{A)^ + 972i?(A)3 + 297R{Af + 36i?(A) + § 

32i?(A)3 + 1 

m{R,A) = 

54/2(A)4 + lAIAR{Af + ^R{Af + 54i?(A) - ^ 

32i?(A)3 + 1 

n{R,A) = 

^fR{AY + 6349i?(A)3 + + 22977 

82R{Af + 1 


The solution a(A) has four possible branches. As in the case of the two-matrix model, we 
require the relation, R{z,A) — 1 ~ (A — ^ which determines a(A) uniquely. By 

using a(A), the other functions are given as 


6(A) 

5(A) 

d(A) 

^(A) 

h(A) 

A 


i?(A)^A {1367 - i2h(A) - A?>2a{kf - 64a(A)3 - (144 + 64h(A)) i?(A)} 
32 {2 + 9i?(A) +4a(A)i?(A)} ^ 

Q 97 S'! 

-a(Kf + ^S(A) + 3^1 fl(A) - b{h)R(hfl\ 


R{Af/^ 


2a(A) - 


5(A) = 


R{Af 


/(A) = -i?(A), 


|2a(A) + 5}fl(A)>y 


■S(A)■^-la(A) + ^-2KA)fl(A)‘''^ 


- {2a(A)6(A) + l6(A)| fl(A)‘/2 + |2a(A) + H| fl(A) 
+ |3a(A)" + ya(A) + aLl R{h_f + 2i(A)R{Af/\ 


(A.17) 


The classical functions Xi(^,A) and X 2 {z,A) enable us to evaluate the expectation 
values -^{trA^) (= -^{trC^)) and For example, can be determined by 


A 

TV 


{B-) 



It is hard to calculate the expectation values ^{tr{A^B^)) and ^{tr^A^B^C^)) directly. 
These expectation values, however, can be reduced to -^{trA^) by using some kinds of 
the Schwinger-Dyson equations. 
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